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Abstract 


The Kerman-Klein-Donau-Prauendorf (KKDF) model is a linearized ver¬ 
sion of the non-linear Kerman-Klein (equations of motion) formulation of the 
nuclear many-body problem. In practice, it is a generalization of the stan¬ 
dard core-particle coupling model that, like the latter, provides a description 
of the spectroscopy of odd nuclei in terms of the corresponding properties of 
neighboring even nuclei and of single-particle properties, that are the input 
parameters of the model. A divers sample of recent applications attest to 
the usefulness of the model. In this paper, we first present a concise general 
review of the fundamental equations and properties of the KKDF model. We 
then derive a corresponding formalism for odd-odd nuclei with proton-neutron 
number (Z, N) that relates their properties to those of the four neighboring 
even nuclei {Z -\-1, N -\- 1), {Z — 1, N -\- 1), {Z -\-1, N — 1), and {Z — 1, N — 1), 
all of which are required if one is to include both multipole and pairing forces. 
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We treat these equations in two ways. In the first, we make essential use 
of the solutions of the neighboring odd nucleus problem, as obtained by the 
KKDF method. In the second, we relate the properties of the odd-odd nu¬ 
cleus directly to those of the even nuclei. For both choices, we derive equa¬ 
tions of motion, normalization conditions, and an expression for transition 
amplitudes. We also resolve the problem of choosing the subspace of physical 
solutions that arises in an equations of motion approach that includes pairing 
interactions. 
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I. INTRODUCTION 


The Kerman-Klein-Donau-Frauendorf (KKDF) model for odd nuclei was introduced and 
applied [1-6] as a semi-phenomenological approximation to the Kerman-Klein (KK) self- 
consistent formulation of the equation of motion approach to nuclear collective motion [7-11]. 
As such it generalizes phenomenological core-particle coupling models, to which it can be 
shown to reduce in various limits [12], The past decade has witnessed further development of 
the theory and additional applications [13-21] including, for example, a suggested solution of 
the Coriolis attenuation problem [17,18]. A review of this more recent work is in preparation 
| 22 ], 

The main purpose of this paper is to show that a formalism of the KKDF type, at the 
same level of completeness as for odd nuclei can be constructed for odd-odd nuclei. A hrst 
important step in this direction has already been made by Starosta et al who have applied a 
restricted version of the formalism to the phenomenon of chirality in odd-odd triaxial nuclei 
[23]. The restriction is the ommission of pairing interactions. When the latter are included, 
we face, among other difficulties, the problem that the manifold of solutions is four times 
the size of the manifold of physical solutions. More recently Koike et al [24], these authors 
have applied an approximate form of the formalism developed in Sec. III. 

As a preliminary step, in Sec. II, we review the KKDF program for odd nuclei. We 
do this in a form which is both more general and more concise than can be found in our 
previously published work, and which sets the stage for the work on odd-odd nuclei that 
follows. It is more general in the sense that the equations are not restricted to deformed 
nuclei. It is more concise in the sense that in our published work, we have described up 
to three different methods for choosing the physical subspace of solutions, whereas here we 
choose that one of these methods that should work in all cases and is, in any event, the 
simplest to implement. 

In Sec. Ill, we present the first of two methods that can be used for odd-odd nuclei. 
We refer to this as the sequential method in that it solves the problem by two successive 
applications of the KKDF approach to odd nuclei, utilizing the solutions for neighboring odd 
nuclei to derive equations for an odd-odd nucleus relative to it neighboring odd nuclei, so 
that the method involves only single-particle coefficients of fractional parentage (CFP). In 
Sec. IV, in an approach that treats the pair of odd particles symmetrically, we derive a set of 
eigenvalue equations and attendant orthonormalization conditions for two-particle (proton- 
neutron) coefficients of fractional parentage. These amplitudes relate the given odd-odd 
nucleus to any of four neighboring even nuclei. For both approaches, we solve the problem 
of choosing the physical subspace of solutions. Finally we derive for each case formulas for 
single-particle transition matrix elements that clearly separate collective and single-particle 
contributions. 
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II. REVIEW OF MODEL FOR ODD NUCLEI 


A. Equations of motion 

In this section we shall derive a version of the Kerman-Klein (KK) equations based on 
the Hamiltonian (2.1) given below. These equations, when taken literally, define a non¬ 
linear problem for the self-consistent study of the properties of an odd nucleus and of its 
immediate even neighbors. However, the version of the theory developed here, referred to 
as the Kerman-Klein-Donau-Frauendorf (KKDF) model, has a more modest goal. This goal 
is achieved by making such further approximations as to reduce the problem to a linear 
eigenvalue problem for the properties of odd nuclei, assuming the required properties of the 
neighboring even nuclei to be known. This can be done only if the Hamiltonian can be 
chosen of sufficiently simple form that the matrix elements of its ingredient multipole and 
pairing operators can be related to observed properties of the even neighbors. Even with 
such simplification, the resulting theory generalizes previous core-particle coupling models. 

We start with a shell-model Hamiltonian of the form 

1 1 
a ^ 

abed LMl 

+ 9 X! X! G'a6cd(^)^LMi(«^)"4LMi:,(cd). (2.1) 

^ abed LMl 

Here ha are the spherical single-particle energies referred to the nearest closed shell, a refers 
to the standard set of single-particle quantum numbers, including in particular the pair 
{ja, ^a) and a refers to the same set with rua omitted. The charge conservation requirement 
means that only the matrix elements of interactions F and G which fulfill the condition 

qa + qb = qc + qd, ( 2 . 2 ) 

where qa is the electric charge of a nucleon with the set of quantum numbers a, do not vanish 
and enter in the Hamiltonian of Eq. (2.1). In the KKDF model we assume additionally two 
more restrictive conditions for the interaction matrix elements, namely: 

1. the charge exchange interactions are excluded, i.e. qa = qc and qb = qd for non 
vanishing matrix elements Fa^ysp, 

2. only the pairs of like nucleons are correlated, i.e. qa = qb and qc = qd for non vanishing 
matrix elements Ga/S'yS- 

B^lml particle-hole multipole operator, 

Sfsijamajb-mb\LML)ala0 

rriamb 

= (2.3) 
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and ^lml particle-particle multipole operator, 

H {jamajbmb\LML)aial, (2.4) 

rriamt, 

where {jimij 2 m 2 \j'm) is a Clebsch-Gordon (CG) coefficient, = (—1)'^““™“. The coeffi¬ 
cients F are the particle-hole matrix elements. 


Facdb{L) = ^ S^Sp{jamajc “ mc\LML) 
m's 

X {jdrudjb - mb\LML)Fa^sf 3 , (2.5) 

and G the particle-particle matrix elements 


Gabcd{L) = Y^{jamajb'mb\LML) 
m's 


X {jcrricjdmd \ LM^Gap^s ■ 

(2.6) 

Assuming the matrices F and G are real, we have 


Facdbi^L') Fdbaci^Ly 

(2.7) 

Gacdb{L) = Gdbac{L) 




= i-iy^^^‘^-^+^Gacbd. 

(2.8) 


The task is to obtain equations for the states and energies of an odd nucleus assuming that 
properties of immediately neighboring even nuclei are known. The states of the odd nucleus 
(particle number A) are designated below as \J^v), where v denotes all quantum numbers 
besides the angular momentum J and its projection p. The states of the neighboring even 
nuclei with particle numbers (Ail) are written, in a parallel notation, as |/M?7.(Ail)). The 
corresponding eigenvalues are Ejy and respectively. We hrst obtain the operator 

equations of motion (EOM), bar indicating reversal of the sign of the single-particle magnetic 
quantum number. 


[Uq., hf] h^Cla 

+ EE 

bd-i LM 

+ EE 

bd-f LM 

[al, H] = -h”al 

-EE^ Ma'ITT'ajc “ 771^ LM) a))Facd6(A) (d6) 

bd'y LM 

- E '^Uarriajc - mc\LM)a;yGacbd{L)A[jy^{bd). (2.10) 

bd'y LM 


Here 
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Kcdb = ^iFacdb + = Facdb, ( 2 . 11 ) 

hg = hg — FcacajL) , (2.12) 

^ Ljc + J- 

K = ha + Y. |^^(2Gaeac + (2.13) 

In consequence of (2.11), we may replace F by F . 

The appearance of different single-particle energies in the two equations may be traced 
to the rearrangement of operators required to have the EOM in a form necessary to achieve 
our aims. This requires, as we shall see below, that the multipole and pairing operators 
occur on the extreme right. The matrix elements of these equations provide expressions 
that determine the single-particle coefficients of fractional parentage (CFP), 

Vj^u{c(IMn) = {J fj,i'\aa\IMn{A + 1)), (2.14) 

U j^y{aIMn) = {J ^v\a\\I Mn{A — 1)). (2.15) 

To hnd equations for these quantities, we form the necessary matrix elements of the EOM 
and evaluate the interaction terms by inserting the completeness relation for the states of 
the appropriate even nuclei between the single-fermion operators and the multipole or pair 
operators. 

In terms of a convenient and physically meaningful set of energy differences and sets 
of multipole helds and pairing helds dehned below, we thereby obtain generalized matrix 
equations of the Hartree-Bogoliubov form 


SjyVj^u{aIMn) 

= (d + 

F^aIMn,^I'I FI Tl ), 

£juUj^u{aIMn) 

= { — M'n') 

M Tl). 

Here 

Sj, = -Ej, + 

(K-Xa), 

- A'A, 

(Aibl) _ XXX X d 

^alMn.'yl'M'n' — (Ja-fOu'CMM^Onn'K-^ln ~ -^0 

^ al Mn^'jF M'n' EE^ ^{.janiajc - rriclLML) 

L bd 

xFagdb{L){l'M'n\A ± l)|HiMjd5)|/Mn(H ± 1)) 

^a.IMn^'yl' M'n' E E Uarriajcmc \ LMl) 

L bd 

xGgML){FM'n\A-l)\ALMAdb)\IMniA + l)). 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 
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Furthermore refer to the ground state energies of the neighboring even nuclei, the 

matrix elements of F^ are derived from those of (2.22) simply by the replacement of the 

operator B hy B\ and the matrix elements of Al are similarly derived from those of A by 

the replacement of A by together with the interchange A ± 1 ^ A ^ 1. Finally is 
obtained from by the replacement of by h^. 

To specify fully solutions of the equations given above, we must develop orthonormal¬ 
ization conditions for the CFP that fix their scale. Orthogonality conditions can be derived 
from the equations of motion themselves. A normalization condition, on the other hand, is 
obtained by taking a suitable matrix element of the summed anticommutator, 

^{a„,at,} = 0, (2.24) 

Q. 

0 = ^(2j, + l). (2.25) 

3 a 


We thus find 

^ E + \Vj,M-IMn)\^] = 1. (2.26) 

aIMn 


B. Equations for reduced matrix elements 

To apply the Wigner-Eckart theorem to obtain the EOM for the reduced matrix elements, 
we utilize the following dehnitions for the latter (which suppress nucleon number): 


Vj^MlMn) = (2.27) 

Uj^y{aIMn) = {IMjama\J^i)ujy{aIn), (2.28) 

{l'M'n'\BLMAbb')\IMn) = (-l)^-^^(/ML - Ml\I'M'){ l'ri\\BL{bb')\\In), (2.29) 

{I'M'n'\ALMAbb')\IMn) = - Ml\I'M'){I' n'\\AL{hh')\\In), (2.30) 

{I'M'n'\B{j^^{bb')\IMn) = {I M LMl\I' M'){I'n'\\B{{bb')\\In), (2.31) 

{I'M'n'\A{j^^{bb')\IMn) = {IMLMl\I'M'){I' n'\\A{{bb')\\In). (2.32) 

Assuming the reality of the multipole and pairing matrix elements, we also have 

{I'M'n'\BLMAlMn) = {!'M'LML\IM){In\\Bl{bb')\\I'n'), (2.33) 

{I'M'n'\ALMAlMn) = M' LML\IM){In\\A{{bb')\\I'n). (2.34) 


With the help of these dehnitions, we can transform Eqs. (2.16) and (2.17) into the forms 

£juVju{aIn) = + 

+ ^ I'n')vjv{a'I'n') 

a'I'n' 

+ E Aj(a/n|a'/'n')'Uj,y(a'/'n'), (2.35) 

a' I'n' 
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where 


EjyUju{aIn) = (-e" +^'')uju{aln) 

— ^ ^\aln\a'l'n')uj,y{a'I'n') 

a'I'n' 

+ ^ A^j{aln\a l'n')vji,{a'I'n), 

a'I'n' 




A j{aln\a' I' n') 


vY-^\aIn\a'I'n') 


/^\{aln\a' I' n') 


^,-iy.+J+/^(2/,+ i)(2L+i) 

Lbb' I 

xK,,fefe,(L)(/V||Si(b6')||/n), 

^(_l)ia+J+/+y(2j/ + i)(2L + 1) 

Lbb' 

xGww(i)(/n|l4(i>'>')IUV). 

J-(_l)i.+x+/^(2;,+ i)(2l,+ i) I 

Lfefe' I 

xFaa'bb'{L){I'n'\\Bl{bb')\\In), 

^(_l)ia+J+/+y(2j/ + i)(2L + 1) 

Lbb' 

xGaa'bb'{L){I'n'\\Ai{bb')\\In). 


I r L) 

ja' ja J J 



I r L) 

ja' ja J j 



The normalization condition (2.26) becomes 

^[|nj, 2 (a/n)|^ + \uj,y{aln)\‘^] = 12. 

aln 


(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 


The equations derived above dehne a linear eigenvalue problem, provided we supply from 
the outside the single-particle energies ha, the reduced matrix elements of the included mul¬ 
tipole and pairing forces, and the excitation energies of the neighboring even nuclei. In the 
underlying (self-consistent) theory these quantities, other than the single-particle energies, 
can themselves be expressed in terms of the CFP v and u. In practice, characteristics of 
even nuclei expressed in terms of the reduced matrix elements of single-particle operators 


FlMl — fac{L)BLM,,{<^c) 

ac 

(2.42) 

and pair transfer operators 


Glml = '^lab{.L)A\^^{ah) 

(2.43) 


ab 


are available rather than the reduced matrix elements of two-body interactions. To make 
use of them in the equations (2.35) and (2.36) it is convenient to present the interactions 
appearing in Eqs. (2.37) - (2.40) as a sum of separable interactions of the form; 


Facdb{L) — —KL{(laQb)fac{L)fclb{L), (2.44) 

GabcdiL) = -gL{qa)lab{L)'^cd{L). (2.45) 


Then the interactions are parameterized by a few strengths k,l and ql which can be either 
htted to the experimental data or estimated theoretically. 



C. Physical solutions 


The equations that we have derived have the form of generalized Hartree-Bogoliubov 
(HB) equations. We summarize the content of Eqs. (2.35) and (2.36) in the condensed form 

•=(:)■ 

_ / e' + A 

\ At -e” + - rt(^-b 

The HB strnctnre of these eqnations implies that only half of the solntions refer to 
physical states. In the standard gronnd-state problem, the solntions divide into two sets with 
reversed energies, the positive energies representing the physical solntions. The solntions of 
Eq. (2.46) do not divide so neatly. The resolntion of this dilemma starts by identifying a 
piece of the Hamiltonian Ti that has snch a simple property and then initially to ”tnrn off’ 
the remainder of the operator. This is done with the aid of the orthogonal matrix C that 
interchanges particles and holes. 


(2.46) 

(2.47) 

(2.48) 


and its transpose C, and by dehning the operator 

= ^(77 - cnc) 

{no)n i(A + At)\ 
i(A + At) -(7f.)n 

(7fo)ii = ^(e' + e") + ^(rA+b + rtA-b) + i(a;A+i) _ 



(2.49) 


(2.50) 

(2.51) 


Because 


CHoC = -Ho, (2.52) 

if T is an eigenstate of Tio with eigenvalne S, then (74/ is an eigenstate with eigenvalne —S. 
As in the simple case, the solntions with positive eigenvalnes are the physical solntions for 
our limiting case. 

Next we tnrn on the remainder of the Hamiltonian, namely the even part 

n^ = ^{n + cnc), ( 2 . 53 ) 

onr aim being to keep track of the physical solntions. In the applications carried ont to 
date, we have described several methods for carrying ont this program. Initially we de¬ 
scribed methods based on tnrning on the ’’pertnrbation” slowly and following the physical 
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solutions by continuity arguments, but in the end realized [16] that by invoking the no¬ 
crossing theorem for two levels of the same angular momentum, we can simply identify that 
half of the solutions of a given angular momentum with the largest energies as the physical 
solutions. From a practical point of view, it suffices to diagonalize the Hamiltonian Tig 
using the complete set of states {physical and spurious) generated by Ho and selecting the 
largest half of the eigenvalues as the physical solutions. The diagonalizationof He within 
the subspace of physical (positive energy) states of Ho performed originally when solving 
the model [3] can lead to a bad approximation of physical solutions of H or even give some 
unphysical solutions, since matrix elements of He between physical and unphysical solutions 
need not be small. 

It is of interest to contrast this procedure with the one used earlier in which the term 
He was turned on adiabatically and the physical solutions followed by using a wave func¬ 
tion overlap argument. This procedure is based on the assumption that the physical wave 
functions change slowly during such a procedure. It is precisely this assumption that fails 
in the neighborhood of an avoided crossing, because when this occurs, it is well-known that 
there is an interchange of wave functions between the two levels involved. In other words, as 
opposed to the simple argument based on the ordering of the levels, the set of wave functions 
assigned as physical must be modihed as one passes a near crossing. 

This brings us to another issue that is both technical and physical. The simplest ap¬ 
plication of the KKDF method is to cases where there is well established band structure, 
either rotational or vibrational, of the same type for both neighboring even nuclei. The 
problem is then to classify the states of the odd nucleus into bands. For this case, the study 
initially of Ho can be useful. This is because for the states belonging to the same band, 
states of different J are practically degenerate, because of the smallness of 
This was the method used in our early work [13-16]. For more complicated situations, we 
can identify different band members by the structure of the states, in the sense that the 
expansion coefficients in terms of a given basis of states vary slowly with angular momen¬ 
tum [20]. Consistent with the identihcation by state vector, we should equally be able to 
associate states into bands by calculating transition rates of a suitable collective operator, 
usually the electric quadrupole operator. 


D. Matrix elements of single-particle transition operators 

We complete the exposition of the general formalism for present purposes by deriving 
formulas for transition amplitudes of a general (charge-conserving) one-body operator. We 
choose this operator to be a tensor of rank L, Tlm^, that we write in the form 

Tlml (2-54) 

The notation is such that the quantities tap include a product of matrix elements of single¬ 
particle operators and of associated coupling strengths (charges, gyromagnetic ratios, etc.) 
We wish to calculate the matrix element {J'p' v'ITlmlIJ■ To carry through the cal¬ 
culation, we substitute for the ket a formally exact expression in terms of the action of 
single-particle operators on the states of the core [25,5,6] 
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|jyuz/) = ^ Y. [U.TrJaIMn)al\ IMn) 

aIMn 

+yjfMiy{oiIMn)aa\IMn)], (2.55) 


where an underline identifies the lighter of the two cores and an overline the heavier one. 
That Eq. (2.55) represents an orthonormal set can be proved by first showing that the 
orthogonality of different states follows from the equations of motion (2.16) and (2.17) and 
then showing that the normalization follows from the CFP normalization condition (2.26). 
By using the commutation relations and completeness, this leads to the following expression 
for the transition element: 

alMnl'M'n' 

x( I'M'n' \TT.M,\IMn) 

+Vj,^,,iaI'M'n')Vj^,{aIMn){TW\TLM,\^)] 

+ Q E/ {al Mn)Ujfj_u {a'I Mn) 

a^a\IMn 

-Vj^^{aIMn)Vj>^iyi{6iIMn)]. (2.56) 


This is now evaluated by use of the Wigner-Eckart theorem with the help of Eqs. (2.27) 
and (2.28) and the following additional definitions of reduced matrix elements: 


x{J'u'\\Tl\\Jv), 

iiy-M 

{I'M'n'\TLM,\IMn) = ^^E_ (Tm7 - M\LMl) 
x(JV||T^||/n), 

^ — lyc-rric 

ta'y y ‘2L I 1 'rflc\LM 


(2.57) 


(2.58) 

(2.59) 


We thus find 


x-^^/(2J+ 1)(2J' + l)[n7,Ta/n)M7',/(a/V)(/V||Tr||/n) 

+vj^ {al n) vj>yi {al'n') ( I'n' 11 Tf, 11 /n)] 

+ ^ E taa'\l{2J + 1)(2J' + l)[(-l)^’“'+^+-^'+^ I ^Auj,,yaln)ujya'ln) 

aa'In \ ) 

+ (_l)Ja'+^+<^'+i I ^ I vjyaln)vji^ya'ln)]. (2.60) 


We thus have a clear separation into collective and single-particle contributions. 
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III. ODD-ODD NUCLEI 


A. Equations of motion 


We turn to the problem of deriving a general core-particle coupling model for odd-odd 
nuclei analogous to the model derived for odd nuclei in Sec. II. Given an odd-odd nucleus 
with Z protons and iV neutrons, we shall relate its properties to those of four neighboring 
even nuclei with proton-neutron numbers {Z + 1,N +1), {Z +1,N — 1), {Z —1,N +1) and 
(Z — 1, A — 1) respectively. In the following development, we shall continue to use Greek 
letters for a general single-particle level, but shall use p, p', etc. to indicate proton levels 
and n, n', etc. to indicate neutron levels. To relate the properties of the target odd nucleus 
to its four even neighbors, we need the equations of motion for four pairs of operators that 
we present for initial convenience in an uncoupled form. 


Pp(^ni hf] 




ClpQin, Hi 



~^F'pp'ijp'OjpiCLni^dpiOjfj) -\- ^prClpCl^' 

~\~Gpp'p"p"'ClpiClniflp'"0,p"^ T Gnn'n"n"'OipO'ri'^^iT-'" 

^pp'n'n^p' ^n' i 
[hp h^(lpCl\^ 

FFpp!f)fji(lpio)^(y(lpi(lij) 

FGpp'p'Ip'll(XpfO}^{jXpf"Clpf'^ Gnnin"n'"^p^n' 

T FppirfiYi' (^p' > 

( hp T h^'jOipCifi 

^pp'PIS'^p'^nifl T Ffifi'Q^pdfii {o,pidpl 

Gpp'p"p'"Clpidfii^Clp/iClpiii^ T GfiYi'n"n"'(Sj^p(S^i(^dYi'" 

FFpp 

'nn' i 

i-hp - K)alal 

^pp'p''p'" ^p' Gnn'n"n'"^p0^n'{ 

Fpp'n'n^pf^n' ? 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


Notice that we have not inclnded nentron-proton pairing interactions. 

We shall stndy matrix elements of these eqnations between the states {JMjs\ on the 
left, an inclnded state of the odd-odd nucleus (Z,N) and the appropriate one of the states 
\aTRMpir) of the even nucleus (Zal, Arl), cr = ±, r = ±. For the further development of 
the formalism, in particular the reduction to equations for reduced matrix elements, there 
are, however, several choices. In this section, we shall develop a method that makes maximal 
use of the formalism for odd nuclei, and, as a consequence involves only single-particle GFP. 
This method, which treats the nentron and proton asymmetrically, we shall refer to as 
seqnential coupling. In the next section, we shall develop an alternative, referred to as 
symmetrical conpling, that bypasses any use of the results for odd nuclei. 
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B. Equations of motion in sequential coupling 


By introducing a complete set of states of the appropriate odd nucleus between the 
neutron and proton single-particle operators, we write 

{JMjs\apan\ + +IMjr) = 

= ^ XjMjs{pJnMnrn)v}^lf^^^{nIMir), (3.5) 

Jn I^n ^ n 

{JMjs\apal\ + -IMjr) = Mir) 

= ^JMjs{pJnMnrn)UfjM^^^{nIMir), (3.6) 

Jn Mji V n 

{JMjs\alan\ - +IMir) = 

= H YjMjsipJnMnrn)v}~li^^^{nIMir), (3.7) 

Jn A/ri V n 

{JMjs\al,ai\ - -IM,r) = >I'|,y’,(pn/Af,r) 

= E YjMApJnMnrn)U<,-i,^,JNIM,r), (3.8) 

Jn^n'f'n 

Here 

= {aJnMnrn\an\(7 + JM/r), (3.9) 

= {aJnMnrn\an\(7 - IMjr), (3.10) 

are two sets of CFP amplitudes for odd neutron nuclei, which can be calculated using the 
formalism for odd nuclei developed in Sec. II. On the other hand the amplitudes 

XjMjsipJnMnVn) = {JMjs\ap\ + (3.11) 

YjMjsipJnMnVn) = (JMjs|aj,| - JnMnVn), (3.12) 

are single-particle CFP relating odd and odd-odd nuclei. The aim of the present coupling 
scheme is to obtain equations to determine the amplitudes X and Y. Before proceeding along 
these lines, we remark that there is a related sequential scheme obtained by starting with 
two-particle amplitudes in which the order of the single-particle operators is interchanged. 

The next step is to write out equations for the amplitudes Mjr). To fix 

additional notation, which will be understood immediately to be a modihed form of the 
notation of Sec. II, we exhibit just one of these equations (using the summation convention), 

{-Ejs + E++ir)¥fMl{pnIMir) = {h^ + h'j¥/j^l{pnIMir) 

+T^++\pjMir\p'I'Mi,r')'^^^jy^l{p'nI'Miy) + A^+\pjMir\p'I'Mi,r')'^^j~Ml{p'nEMi,r') 
+TY+\nJMir\n'EMrr')'^^fMl{pn'I'Mi,r') + A^+\nJMir\n'I'Mi,r')'^^fMl{pn'I'MiM) 
-Epp>n>n'^^:^Mjs{pnIMir). (3.13) 

Into this equation and its three partners, we substitute Eqs. (3.5)-(3.8) and recognize that 
the result can be simplihed by the use of equations such as 
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= r^++\nIMir\n'EMi,r')vj^l,^,^{n'E 

+A(+) {nIMir \ n'l'Mp r') Uj^lf-nrr. {n'I'Mr r') (3.14) 

and its partners. The resulting equation for is then combined with the corresponding 

equation for , contracting the first with a factor and the second with a factor 

to as to permit use of the normalization condition (2.26). We carry through a corresponding 
procedure for the pair of amplitudes and 

We thus obtain the pair of equations 

( Ejg -f- E^j^^^^Xjj^jgi^pJnM^Tn) = hpXjMjsi^pJnMjiTri) 

+r^ ^ (^PJ^M^TyiIp Jn'JMjsip' Jn'n'^ 

-\-X{j)JnMYiTn\p Jn' MniTn'yEjMjgip Jn' M^'Tn'^ 

+ V^"^\pJnMnrn\p'Jn'Mn'rn')XjMjs{p'Jn'Mn'rn'), (3.15) 

( Ejg + E_j^^^yYj]r[jg{j)JnM^T^ hpYjjnjgijjJnMjyVn) 

ph ^ [j)Jj^M^Tji\p Jn'^n'f'n'y^JMjsip Jn'^n'^n'^ 
(pt/„iVf„r^I p Jji'M^'T yi'^XJ^ jg{j) .J^iM^iTn'^ 

+ V^~\pJnMnrn\p Jn' Mn'rn')YjMjs{p' Jn' Mn'Vn'). (3.16) 


Here 


nlMirl'Mj/r' 

><[bi.rJ"«f/0r‘++>(P^Wr|p7'M,,r')q;P,,_,(n/'A/,,r') 


^{^pJnMnfnlp Jn'Mji'Tfii^ — ^ ) 

nlMirl'Mj/r' 

><[b,Tr„("«f,r)Al+)(p/M,r|p7'A/,.r')q;l4,.,,(i7'M,,r') 

V’^Yi^pl^MnrnW Jn'Mn'rn') = 

r'^^“^(pJnM„r„|p'Jn'M„/rn') = Y 7^ 

nlMjrl’Mj/r' 

><lb:hr.(>>M*>-)r*'-+'(p/M,r|p7'Afpr')q;L.,r.,(n/'Afrr') 

1 

X\pJnMnrn\p Jn'Mn'Vn') = Y 


nIMjrPMjir 


I ^(n) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 
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V^~\pInMnrn\p' Jn'Mn'Vn') = 7^ 

X l^tLrr. Mjr) Fpp^nn' 

(n7M7r)Fpp/,,/]. 


(3.21) 


(3.22) 


The correctness of the above equations can be verihed independently. By starting with 
Eqs. (2.9) and (2.10), we can derive equations of the form (2.16) and (2.17) with potentials T 
and A that refer appropriately to the odd systems, rather than the even system and with no 
overt sign of the neutron-proton interaction terms. These equations are readily transformed 
into the results given above by the application of Eq. (2.55), just as the latter was applied 
in Sec. II.D to express transition matrix elements between odd states in terms of matrix 
elements between even states and single-particle CEP. 

The hnal goal of this section is to obtain equations of motion for the reduced matrix 
elements. For this purpose the only dehnitions needed to supplement Eqs. (2.27)-(2.32) are 
those for the reduced CEP relating the odd to the odd-odd nuclei. 


XjMjsiPJnMnTn) = (“l)^'’ {JnMnjpmp\J Mj)xj sijpJnTn) , (3.23) 

YjMjsiPJnMnTn) = ( j^TTlp| JMj)?7j^( jp J„r„) . (3.24) 


It is then a straightforward exercise in angular momentum algebra to derive the equations 

{~Ejs + E^j^rn)XJs{jpJn'^n) — hpXJs{jpJn'^n) 

^ i^jpJn'^ n\jp'Jn''^n')Xjs i^jp' ) 

^ Op J'nTn\ jp' VJs (,jp' 

^ {^jpJn'^n\jp' Jn^'^n'^Xjs{^jp^ •) (3.25) 

( Ejg + E—j^^^^Tljg(^jpJnTn) hpTjjs^jpJn'^n) 

^^nO VJs i^jp' 

{,jpJn^n\jp' Jn'^n'^Xj s i^jp' Jn'^n ') 

^ {^jpJn'^n\jp'si^jp' ^ (3.26) 


where 


r^ ^ (Jp Jn'^n\jp' 


^(Jp Jn'^nl jp' 


^ \ Jn jp \ ^ \ 

^(n) \ jp' J \ ^ ^ } 

x(-l‘)7+J.+P+''+i-i-.y(2L + 1)(2J„ + 1)(2J„. + 1)(2/' + 1) 

3 j Jn jp \ ^ jn 1 

fl(n) \ jp' '^n' J \ Ti' I ^ J 
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(3.28) 


^ ^jpJn'l'n\jp' Jn'f'n'^ — 


x(_i)ip+i"+-^+^'-i^(2L + 1)(2J„ + 1)(2J„, + 1)(2J' + 1) 
X Gpp^aa'iL) [v^Zl„ {jJr)v^j~Jr^, {jJ'r') (/V'| \A]^\aa') \\Ir) 
+«£n ijnlr)u^j~jr^, {jni'r') {tta') 11 Jr)], 

J jp 1 / jn' jn L 1 


^ i^jpJnf’nlip' 


f^(n) \ jp' ^ j Jn' ^ j 
x(-i)>+i-+'^-+'^+^(2L + 1)^(2J„ + 1)(2J„/ + 1) 
^ ['’^JnTn n'I v) Fppt n'n{,L) 

+ ( -1 WjtL {3nlr)u^^]r^, {jn' I r) Fpp,nn'{L) ], 

1 / Jn jp 1 / ^ Jn jn 1 


(3.29) 


^(n) \ jp' J'u' J \ ^ ^ ) 

X (_i‘yp+Jn+J+/'+i-L^(2L + 1)(2J„ + 1)(2J„, + 1)(2J' + 1) 
X Fpp'aa' [vtl i3nIr)vFX, {jj'F) ( JV'11 11 

+4tl(jn/r)4+) ,(j„JV)(/V||J?i+-Va')l|/r)], 


(3.30) 


i^jpJn^nljp'Jn'^n') 


^ j Jn jp 1 / ^ Jn jn 


^ i,jpJn^n\jp' Jn'^n'} 


f^(n) \ jp' Jn' L j y Jn' I L j 

X {-iyv+jn+J+i'-L^^2L + 1)(2J„ + 1)(2J,, + 1)(2J' + 1) 
X Gpp^aa' ijnM'>^jll„, UnU') (JV'11 {att') \\Ir) 

+41 (Jn/r)4+t^, (j.JVO(/VI |4^-VaO I I/O], 

^ j Jn jp \ / jn' jn L 1 


(3.31) 


f^(n) \ jp' J'u' J j Jn' I J 
x(-i)ip+in'+Jn+J+i(2L + 1)^(2J„ + 1)(2J„, + 1) 
X [ 41 (in/r)H (in'/r)Fppw(L) 

+ (-1 y ■■^“1 Jin/r) 4“1 (j„, Jr) Fpp/„,„ (L)]. 

The normalization condition associated with this formalism is 


(3.32) 


^ ' [IXjs(jp'Ai'^n) I T l'JJs(jp'Ai'^n) I ] ^(p)- 

jpj—nrn 


(3.33) 


C. Physical solutions 

The problem of choosing the physical solutions of Eqs. (3.25) and (3.26) can solved by 
simply repeating the arguments given in Sec. II.C. This is seen immediately if we rearrange 
the energies in these equations so that they resemble exactly the corresponding equations 
(2.35) and (2.36). We thus write 
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^JsXJs{jpJn^n) {^p ^+Jnrn)XJs{jpJn'>'n) 


^ (,jpJn^n\jp'Jn''f"n'^X.Js i^jp' Jn''^n' ) 

“1“ ^ i^jp Jn^n\ jp' Jn''^n') VJs {jp' Jn''^n'^ 

“1“^^ ^ {jpJfi'^ n\jp' Jn'^^n'^^J s {jp' Jn''^n '); 

(3.34) 

^Js'3Js{jpJn'^n) ( Cp ^—Jnrn)VJs{jpJn'^n) 

^ {jpJn^ rL\jp'Jn''^n'^^Js {jp' Jn'^n ') 

^ {jpJn'^n\jp'Jn''^n')XJs{jp'Jn''^n') 

“1“^^ ^ {jpJn^ n\jp'Jn'^n'^'3js {jp' Jn'^n') •> 

(3.35) 

^Js = ~Ejs + 

(3.36) 

E = fc - i(£+ - E.), 

(3.37) 

^±Jnrn. -^±J„r„ E±- 

(3.38) 


Here E± are the ground states of the heavier and lighter odd neutron nuclei, respectively. 
With these forms one has an exact parallel to Eqs. (2.35) and (2.36), and thus the arguments 
for choosing physical solutions can be repeated without modification. 


D. Matrix elements of single-particle transition operators 


The result we want can be read off directly from Eq. (2.60) if we replace, appropriately, 
the single-particle CEP v and uhy x and i] and the reduced matrix elements of the transition 
operator between states of even nuclei by the corresponding matrix elements between 
states of the appropriate odd nuclei. We thus obtain 


{J's'\\Tl\\Js) = 


('^) jpJnTnJ^lT^I 




+J„+L I '^n •Jn‘ 


tJ Y) J n ^ 


J' J Jp 


^ '\j n)'^J' s' i^3pJn''^n'^ ^ Jn''^n' \ \ Tl \ \ 

~^X.Js {,jpJn'^n)X.J' s' i^jpJn'^n'^ {Jn''^n'\ |^L | | 

+ 7^ H tpp'3 (2J + 1)(2J' + !)[(-!)I ^ \r]j>s'{jpJnrn)VJsUp'Jn 


+(-l) 


J„/+Jn+J'+1 1 3n jn’ L 

\ J' J Jr, 


XJs ^jpJrJn^Xj' s' ijp' JitJu) • 


(3.39) 


IV. FORMALISM FOR SYMMETRICAL TREATMENT OF ODD NUCLEONS 

A. Equations for reduced matrix elements 

With the aid of the Wigner-Eckart theorem and suitable dehnitions, we proceed to the 
transformation of these raw equations to equations for reduced matrix elements. First we 
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introduce a two-component operator for protons, 

= a = +, 

= 0-1: 0 = ( 4 . 1 ) 

and a corresponding pair, a)), r = ±, for neutrons. We then couple the products, 

OpOl = J2{jpmpjnmn\lm)l3Z{pn), ( 4 . 2 ) 

i 

and define reduced matrix elements bYg{pnlRr), 

{JMjs\l3Y{pn)\aTRMRr) = {RMrlm\JMj)bY{pnlRr). ( 4 . 3 ) 

Here \ JMjs) is a state of the odd-odd nucleus {Z, N) and larRM^r) is a state of the even 
nucleus (Zcrl,A^rl). 

In the following, we also require reduced matrix elements of the multipole and pairing 


operators, dehned as follows, 

{arR'MRYlBLM.mWrRMRr) = - Ml\R'Mr,) 

x{R'r'\\By"\bb')\\Rr), ( 4 . 4 ) 

{aTR’MR,r'\Bij^Ybb')\(TTRMRr) = {RMrLMl\R' MR>){R'r'\\By^^\bb')\\Rr), ( 4 . 5 ) 

{a - R’MR,r'\ALMinn')\a + RMrt-) = - Ml\R!Mr,) 

x{R!r'\\A^Y\nn')\\Rr), ( 4 . 6 ) 

{a + R!MR,r'\A\j^{nn')\o - RMrt) = {RMrLMl\R!MR,){R! r'\\Af\nn')\\Rr): ( 4 . 7 ) 

{-TR!MR,r'\ALM{pp)\ + tRMrt) = {-if-^YRMnL - Ml\R!Mr,) 

x{R!r'\\Ay\pp')\\Rr), ( 4 . 8 ) 

{+TR'MRY\Aij^{pp')\- tRMrt) = {RMRLML\R!MR,){R!r'\\Af\pp')\\Rr). ( 4 . 9 ) 


In the dual equations of motion given below, we also introduce in as close analogy as possible 
with our procedure for the odd-nucleus case, various combinations of energies. The energies 
of the odd states will naturally be specified by i7js, those of the four neighboring even nuclei 
by EfjT-Rr: the ground states of the latter by We then introduce the following differences 


S-Js = 

= -Ejs + ^iE++ + E. 

f- 

+ E_+ 

+ E..): 

(4.10) 

O^cftRt - 

E(jTRr Ear 5 



\iE,, 


(4.11) 

^'n = 

-K- 

- 

E+-) 

— 

-E-): 

(4.12) 


-K- 

\(E^. - 

E-+) 

— 


-E..): 

(4.13) 

< = 

= K- 

\{E,, - 

E+-) 

+ 

\(E.. 

-E..): 

(4.14) 


= K- 

- 

E-+) 

+ 

\(E,. 

-E-): 

(4.15) 


= K + 

- 

E+-) 

+ 

\e__ 

8^ 

-E+Y: 

(4.16) 

e'" = 

p 

+ 

II 

- 

E-+) 

+ 

-{E__ 

8^ 

-E+Y- 

(4.17) 
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With the included matrix elements dehned below, we thereby obtain the following equations 
of motion for the reduced matrix elements, 

Sjsby~g{pnlRr) = ^ H{aTpnlRr\a'T'p'n'l'{p'n'l'R'r'), (4-18) 

a'r'p'n'V Wr' 

with values for the non-vanishing matrix elements of the effective Hamiltonian given by 
the expressions 


'H{+ +pnlRr\ + +p'n'l'R!r') 


— pnlRr\ -I— p'n'l'R!r') 


'H{—\~ pnlRr\ — \-p'n'l'R!r') 


n{--pnlRr\ - -p'n'l'R'r') 


H{+ + pnlRr\ — \-p'n'l'R'r') 
H{+ + pnlRr\ -\ — p'n'l'R!r') 

7i(-|— pnlRr\ - p'n'l'R'r') 

— pnlRr \ + +p'n'l'R'r') 
— \-pnlRr\ + +p'n'l'R'r') 

— \-pnlRr\ - p'n'l'R'r') 

H{ - pnlRr\ -\ — p'n'l'R'r') 

- pnlRr\ — \-p'n'l'R'r') 


(e'p + — UJ++Rr)bp^pi6n,n'Si^li6]i^Ji'6ry 

+r^_^^\plRr\p'l'R'r')6n,n' + r^^^\nlRr\n'l'R'r')Sp,p> 

+V^^^'>{pn\p'n')6ijfSR^R>Sr,r', (4.19) 

i^p ~ ^n~ ‘^+-Rr)bp,p'Sn,n'bl^l'SR^R'6ry 

+r[^~\plRr\p'l'R'r')Sny - Ty~\nlRr\n'l'R'r')Sp^p> 
+V^^~\pn\p'n')6i^i>SR^R'Sr,r', (4.20) 

( —Cp + ^n~ ^-+Rr)bp,p’Sn,n'Si^p6R^R'6ry 

-r\^~^\plRr\p'l'R'r')6n,n' + r^^^\nlRr\n'l'R'r')6p,p> 
+V^~^\pn\p'n')6i^i>SR^R>Sr,r', (4.21) 


{ — ^p — ~ ^ Rr)bp,p'Sn,n'bl^l'S 

-T\^-\plRr\p'l'R'r')Sn,n' - rf 

+V^~~\pn\p'n')6i^i>SR^R'Sr,r', 

A^;\plRr\p'l'R'r')Sny, 

A[j^\nlRr\n'l'R'r')6p^p,, 

A^j~\plRr\p'l'R'r')6n,n', 

A^}^\nlRr\n'l'R'r')Sp^p>, 

AP{plRr\p'l'R'r')Sn,n', 

A[p\nlRr\n' I' R'r')6p^pi, 

A^}~\plRr\p'l'R'r')Sn,n', 

Aj \nlRr\n'l'R'r')Sp^p/. 


R,R'br,r' 

\nlRr\n'l'R'r')6p^pi 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 


The remaining matrix elements 7i(crr| — a — t) vanish. 

The effective interactions that occur in the above equations are 


T^J"\plRr\p'l'R!r') 


T^J"\nlRr\n'l'R!r') 


E( 

Lhh' 


.lyp+P 


,+l+l'+L+R+J 


x^{2R' + 1)(2L + 1)(2/ + 1)(2/' + 1) 
xFpp,w{L){R'r'\\B^r\bb')\\Rr), 

'^(^-lyp+R'+L+R+J 

Lbb' 


I R J 
R' I' L 


in ip I 1 

L I' Jp, j 

(4.31) 
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vY"\plRr\p'l' E!r') 


TY"\nlRr\n'l'B!r') 


A^j\plRr\p'l' R!r') 


/\^j\nlRr\n'l'R!r') 


A^y\plRr\p'l'R!r') 


AY\nlRr\n'l'R'r') = 


V^^^\pn\p' n') = 


■\/(+ )(j)n\p'n') 


X 


xK^,,y{L){R'r'\\B^r\hb')\\RT), 

^ ^ ^_ 'Y^3v^3n-\~L^l+V^R.^J 


Lbb' 


X 


xF^^,,y{L){R'r'\\B^^-^^L-{hh')\\Rr), 

(^_\yv+3n'+^+R+J 


Lbb' 


X 


xF^n'bb'{L){R'r'\\B]^^^\bb')\\Rr), 


Lp"p" 


X 


X Gppyyyi) {RV\\ aP (//') 11 /2r), 

_ Yiyp+F'+L+R+J 


Ln"n'" 


X 


X Gnn'n"n'"{L) {R!v'\ | {n"n"')\\Rr), 

_ y ^ ^_ ^yp-\-jii.-\-L+i+i'+R+J 


Lp"p" 


X 


X GppyyyL) {R'r'W Af^ {p"p"') \\Rr), 
. y^ ^_ y3p+3n'+R+R+-^ 


Ln"n"' 


I R J 
R' V L 


I R J 
R' V L 


I R J 

R' r L 


I R J 
R' I' L 


I R J 
R' V L 


I R J 
R' I' L 


X 


X Gnn'n"n'"{L) (R^F \ \ (n"n'") \\Rr), 

-^(-iy'>+^"'+^(2L + 1) 

L 

{t h' i 1 


I R J 
R' I' L 


= ^{-1)A+F+i+L(^2L + 1) 


jp jn ^ 1 

L I' Jn' J 

(4.32) 


jn jp ^ 1 

L V J,, f 
(4.33) 


jp jn ^ 1 

L I' Jn' j 

(4.34) 


jn jp ^ 1 

L V jp f 
(4.35) 


jp jn ^ 1 

L I' Jn' J 

(4.36) 


jn jp ^ 1 

L V Jp, j 
(4.37) 


jp jn ^ 1 

L I' Jn' f 

(4.38) 


(4.39) 
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(4.40) 


X 


ip jn I 


Fpp'nn' (-^) 5 


jn' ip' L 

V^-+\pn\p'n') = ^(-l)V+^"'+'+^(2L+1) 


X 


jp jn I 


Fpp'nn' (-^)) 


jn' jp' L 

V'^ \pn\p'n') = — 'y^{—l)F'^F+i(^2L + 1) 


X 


jp jn ^ 
jn' jp' L 


Fpp'n'ni^L^. 


(4.41) 


(4.42) 


To the equations of motion, we add a normalization condition that can be derived from 
the anticommutation relation 


Qpllo-n; qI.I (4.43) 

pn 

ii(rt = i:(2ip+l). (4.44) 

p 

Rearranging the order of the operators, taking a diagonal matrix element in the state | JMjs), 
utilizing completeness and the dehnitions (4.2) and (4.3), we obtain the expected result 

|6!r;’(!>"«r)|" = fi,rtf!,„|. (4,45) 

r^rjpjnlRr 


B. Physical solutions 

We expect the space of physical solutions to be only a quarter of the total space of 
solutions. With a little care, we can generalize the method used to identify physical solutions 
for the case of odd nuclei. If we examine the Hamiltonian matrix Ti given by Eqs. (4.19)- 
(4.30), we see that it can be decomposed into a sum 

Ti. = Tip + Tin — oj + Vnp, (4.46) 

describing in turn an odd-neutron nucleus, an odd-proton nucleus, an excitation energy 
matrix, and a neutron-proton interaction energy. We initially turn off the last two terms. 
Next we dehne two four-by-four matrices, Cn and Cp, 

C„=(fy. (4.47) 

(4.48) 

where the underlined entries are each two-by-two matrices and the matrix C_ is the particle- 
hole conjugation matrix dehned in Eq. (2.49). The matrices Cp, Cn commute with each 
other. 
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We then observe that the averages 


Tip — ijpTp + CnTipCn), 
'Hn = 2^Tin + CpTinCp) 


(4.49) 

(4.50) 


each have a structure more symmetrical than their individual terms. Thus, the non-vanishing 
elements of 'Hp are (in a condensed notation) 


Tip{+ + I + +) — Tip (-\— I H—) 




Ti, 


+ 1 —1“) — Tip 


-) 


= -i(E; + £;" + r«-+> + r«-->). 


Tip(+ + 1 —^) — Tip{+ 


Ti, 


2" P 

+ ! + +)= 'Hp 


= -(A(+> + A'-)), 


+ 


= -(A«+' + A«-'), 


whereas Hn has the block-diagonal structure 

'n„ 0 


Hn = 


n 

0 H, 


(4.51) 

(4.52) 

(4.53) 

(4.54) 

(4.55) 


in terms of two-by-two matrices, and 

K. = 




(4.56) 


We infer from their structure that it is the barred matrices that form suitable starting 
points for the antisymmetrization that was the essential step for identifying physical solutions 
for the theory of odd nuclei. We thus dehne the matrices Hpo and 'Hnoi 


Hpo =2^Tip- CpTipCp), 

Tino=-{Tin-CnHnCn). 


(4.57) 

(4.58) 


Notice that Cn commutes with Hpo and Cp commutes with Hno, but the two Hamiltonians 
do not commute with each other. We also have 


CnHnoCn 'Hno: 

Cp'HpoCp 'Hpo: 


(4.59) 

(4.60) 

(4.61) 
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To identify the physical basis we introduce an auxiliary Hamiltonian 


HoH) — Hno + sHpo, 

(4.62) 

which has the following properties 


CpCnHoie)CnCp = -Hois), 

(4.63) 

CpHoie)Cp = Hois). 

(4.64) 


A quarter of the set of eigenvectors of TLoiX) = form the physical basis. From Eq. (4.63) 
we see that for each positive eigenvalue, Ejs^s) > 0 of Hois) which enter in equation 

'Ho{£)^jMjs{£) = ^Js{£)^jMjs{£), (4.65) 

there is a corresponding negative eigenvalue, —Sjs{e), with the associated eigenvector 
CpCn^jMjsi^)- Following the standard reasoning of superconductivity theory we reject 
half of the eigenvectors, those belonging to negative eigenvalues as non-physical. However 
that still leaves too many states. 

From Eq. (4.64) it follows that the Hamiltonians 'Hois) and Hois) have the same set 
of eigenvalues, i. e., Bjsie) = Tjs(—e), with the corresponding eigenvectors ^jMjsi—s) = 
Cp^jMjsH)- This has as a further consequence that each eigenvalue of 7Yo(0), Tjs(O) is, 
apart from the 2J -|- 1-fold magnetic degeneracy, additionally two-fold degenerate with the 
two eigenvectors ‘hjMjs(O) and Cp^jMjsiH- Next we solve Eq. (4.65) for 0 < e < 1 and 
a given J. For e > 0, level ^js(O) splits into a pair of levels, labeled Bjs>iH BjsdHi 
distinguishing the larger from the smaller value. We choose the larger of the two eigenvalues 
as the physical one, recognizing its role as the analogue of the sum of quasiparticle energies. 
We then rely on the no-crossing theorem to maintain the ordering of the physical states as we 
increase the value of e to unity, that appropriate to the Hamiltonian HoiB) = Ho- Since the 
number of positive eigenvalues of the latter are even, it is a consequence of the arguments just 
given that, for a given J, the physical solutions are the odd-numbered positive eigenvalues, 
counting from the largest value. 

The separation of Ho from the original Hamiltonian of Eq. (4.46) is achieved by the 
following decomposition of 


H — 'Ho + He — a; -|- V^p, 


(4.66) 


where 


He = 


^iHp 


Hn) + -^CpCniHp -|- Hn)CnCp 


+ '^HpiHp — Hn)Cp + -^CniHn ~ Hp)Cn 


(4.67) 


The physical eigenvectors of H can further be found using again the methods similar to 
those discussed in Sec. HC for odd nuclei. 
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C. Matrix elements of single particle operators 


Finally, we turn to the problem of deriving a general formula for the transition matrix 
element, {J' Mjis'\Tlml\J^js) of a single-particle operator, T^Ml ^ manner analogous to 
the calculation carried out in Sec. IID . For this purpose, we utilize a formula for the state 
\JMjs) (compare Eq. (2.55)) 

^ ^ '^^jMjs{P^RMRr)alal\aTRMRr), (4.68) 

arpnRMjir 

^^jMjsiP'fT'RMRr) = {JMjs\apal\aTRMRr). (4.69) 

Equations (4.68) and (4.69) describe a set of orthonormal states, as follows from the equa¬ 
tions of motion and the normalization condition (4.45). We thus derive the formula 

{J'ARs'IRm.IJMjs) = —2^ y ¥f^\,(pnRMRr)i%lApnR'Mn'r') 

“(p)“(n) arpnRMurR'Mnir' 

x{aTR'MR,r'\TLMLWrRMRr) 

“(p)“(n) 

Tpp'nRM^r 

+ 0 4 H ^j^MlsipnRMRr)¥j,j)^^^,{pn'RMRr)tn'n 

(p) (n) apnn'RMfir 

-7r4^ E 

(p) (^) Tpp'nRMRV 

xV(-iy^-^^+V-^y 

^ '^j"MlsipnRMRr)¥j,^^^^,{pn'RMRr) 

apnn' RMjir 

(4.70) 

To apply the Wigner-Eckart theorem requires, in addition to obvious adaptations of the 
formulas of Sec. (HD), only the additional formula 

'^^jMjsipnR^R^) = Y^{jpmpjnrnn\lm){RMRlm\JMj)ll-Jp{pnlRr), (4.71) 

i 

which combines the contents of Eqs. (4.2), (4.3) and (4.69). We thus hud 

(JV||Ti||Js> = —2^ (_i)J+«'+‘+y(2J+l)(2,/' + l) 

\L(p)\L(n) 

arlLjpjnRrR'r' 

^ i t } (^^'RX||Ti||(Tri?r) 

X {pnlRr)b^j,'^J {pnlR'r') 

I ^ ( ]^')fl+J'+.v+.?n ^V +1 

^{p)^{n) rll'Ljpj^,jr,Rr 2L + 1 
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x^(2i+l)(2;' + l)(2J+l)(2,/' + l)|^ j ‘J'l 

xtpipb[fj^ {pnlRr)b[jJg} {pnlR'r') 

I ^ / ^\R+J'+.h+j„+i+i' ‘^R' ^ 

^(p)^(n) Tll'Ljpj„,j„Rr 2 L + 1 

xy(2(+l)(2/' + l)(2J+l)(2J' + l)|^ j -J'l 

xtn'nb^js~^ {pnlRr)b^J,~) {pn'lRlr') 

I ^ / ]^')fl+J'+v+.?n ^V + ^ 

^(p)^(n) rll'Ljpj ijnRr 2L + 1 


L J f 
R V I 


XsJ{2l + 1)(2/' + 1)(2J + 1)(2J' + 1) 

X {pnlR'r') 

fi(p)^(n) rll'Ljpj^,j„Rr 2L + 1 

xy(2(+l)(2/' + l)(2J+l)(2J' + l)|^ j -J'l 


Jp Jn ^ 

^ Jp' 


Jn Jp ^ 

^ Jn' 


Jp Jn I 

I L jp! 


jn jp I 
I L jjii 


(4.72) 


Once again we have a clear separation into collective and single-particle contributions. 


V. SUMMARY AND CONCLUSIONS 

A linearized version of the equations of motion approach to the nuclear many-body 
problem, considered as a generalization of traditional core-particle coupling models has 
proved its worth in a number of recent applications to deformed odd nuclei. In this method, 
the basic object studied is a single-particle coefficient of fractional parentage (CFP) relating 
the states of the even nuclei to those of a neighboring odd nucleus. 

In this paper we showed how the same general method can be applied to odd-odd nuclei. 
We started with a review of the formalism for odd nuclei, since it plays an essential role in 
some of the considerations that follow. We then showed that there are three possible formu¬ 
lations for the odd-odd case, two of which we label as sequential and a third as symmetrical, 
terms that characterize the way in which we couple an extra neutron (or neutron hole) and 
an extra proton (or proton hole) to nearby even nuclei, treated as cores. First we study in 
detail the case where we initially couple the odd neutron to the even cores, an example of our 
method for odd nuclei. We then couple the odd proton to the odd neutron nuclei, introduc¬ 
ing new CFP for this relationship, and making essential use of the odd neutron calculations 
for energies and CFP. The second sequential method, not discussed in detail, reverses the 
order of the odd-particle couplings. In the symmetrical coupling, we hrst couple the two 
odd particles together and study directly the relationship of the odd-odd nucleus to the core 
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even nuclei by means of two-particle CFP. In principle all three methods are equivalent, but 
in practice results will differ owing to the need to approximate. In this regard, the existence 
of alternatives that may be compared may be of some practical advantage. 

Because of the presence of pairing interactions the equations for the odd-odd case yield 
four times as many solutions as are physical. In the sequential method, the problem of 
choosing physical solutions can be solved by sequential use of essentially the same method 
as for the odd case. For the symmetrical coupling case, a more elaborate method has been 
devised. 

Concerning applications, approximate versions of the sequential method have already 
been carried out [23,24]. The symmetrical approach remains to be tried. 
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